Inertial effects in nonequilibrium work fluctuations by a path integral approach 
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Inertial effects in ffuctuations of tlie work to sustain a system in a nonequilibrium steady state 
are discussed for a dragged massive Brownian particle model using a path integral approach. We 
calculate the work distribution function in the laboratory and comoving frames and prove the 
asymptotic fluctuation theorem for these works for any initial condition. Important and observable 
differences between the work fluctuations in the two frames appear for finite times and are discussed 
concretely for a nonequilibrium steady state initial condition. We also show that for finite times a 
time oscillatory behavior appears in the work distribution function for masses larger than a nonzero 
critical value. 
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I. INTRODUCTION 

In recent years, fluctuations in nonequilibrium systems 
have drawn considerable attention to a new kind of fluc- 
tuation theorems. These fluctuation theorems are asym- 
metric relations for the distribution functions for work, 
heat, etc., and may be satisfled even far from equilib- 
rium states or for small systems in which the magnitude 
of the fluctuations can be large. These fluctuation the- 
orems have been proved for deterministic thermostated 
systems [J 0, Q as well as for stochastic systems [1, Q , 
and have also been discussed in connection with the 
Onsager-Machlup fluctuation theory {6^. Moreover, ex- 
perimental confirmations for these theorems have been 
obtained JJl, VI, _12J . It has also been shown that the 
fluctuation theorems include the fluctuation-dissipation 
theorem, as well as Onsajger's reciprocal relations, near 
equilibrium states P, 113] . 

In our previous paper based on a generalization 
of the Onsager-Machlup theory for fluctuations around 
equilibrium to those around nonequilibrium steady states 
using a path integral approach, we discussed fluctua- 
tion theorems for a stochastic dynamics described by a 
Langevin equation. For a Brownian particle driven by a 
mechanical force F{xs, s), the Langevin equation for the 
particle position Xg at time s is of the general form 



d Xg dxg r,, \ I A 

ds'^ ds 



(1) 



with the mass m of the particle, the friction coefficient 
a and a random noise Qg. In our previous paper, as a 
nonequilibrium model we considered a dragged Brown- 
ian particle, in which the mechanical force is given by a 
harmonic force F{xs,s) = —k{xs — vs) with the spring 
constant k and the dragging velocity v. Furthermore we 
mainly considered this model under the over-damped as- 
sumption. This assumption can be used for a dynamics 
on a much longer time scale than the inertial character- 
istic time Tm = m/a, and the dynamical equation under 
this assumption is simply given by neglecting the inertial 



term containing the mass in Eq. ([T]), i.e. by 

dxs 1 X 1 > 

= -F{xs,s) + -Cs- 
ds a a 



(2) 



Equation ^ is much simpler than Eq. ([T]), but informa- 
tion of the system on the shorter time scale than t„j is 
lost in Eq. ([2]). It may be noted that Machlup and On- 
sager already developed their fluctuation theory around 
equilibrium not only for the case corresponding to the 
over-damped case Q but also for the inertial case [1]. 
In our previous paper we discussed also a generaliza- 
tion of the Onsager-Machlup theory for nonequilibrium 
steady states including the inertial term Q. However, 
there we treated only one type of fluctuation theorem, 
the so called transient fluctuation theorem 0, which is 
restricted to equilibrium initial conditions. Another fluc- 
tuation theorem, the asymptotic fluctuation theorem [sj , 
which holds for any initial condition (including a nonequi- 
librium steady state [24| ), was not discussed for inertial 
cases in Ref. Q . Different from the transient fluctuation 
theorem, which is correct for all times as a mathemat- 
ical identity 14], the asymptotic fluctuation theorem is 
satisfied in the long time limit only. However, as we will 
discuss in this paper, a variety of interesting inertial ef- 
fects appear for finite times for a nonequilibrium initial 
condition, before the asymptotic fluctuation theorem is 
achieved. Although there are some results for fluctu- 
ation theorems for stochastic systems including inertia 
(isl [l6| . the asymptotic fluctuation theorem with iner- 
tia has not been discussed fully in connection with the 
Onsager-Machlup theory so far. 

The purpose of this paper is therefore to discuss, in the 
context of the Onsager-Machlup path integral approach, 
inertial effects in nonequilibrium steady state work ffuc- 
tuations, including the asymptotic ffuctuation theorem. 
For these discussions we use the Langevin equation ([1]) 
for a dragged Brownian particle without the over-damped 
assumption. The work distribution function is calcu- 
lated explicitly for any initial condition, and its finite 
time properties are investigated. As an important iner- 
tial effect we show a critical value of mass above which 
the work distribution function shows a time-oscillatory 
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behavior. 

The nonequihbrium work used in this paper is based 
on the generahzed Onsager-Machlup theory, as obtained 
in our previous paper [(J]. In that paper we considered 
two kinds of work in two different frames: (A) the work 
Wi done in the laboratory frame (l) and (B) the work 
Wc done in the comoving frame (c) where the average 
velocity of the Brownian particle is zero in a nonequihb- 
rium steady state. A difference between these two works 
is that Wc includes a d'Alembert-like force, which is ab- 
sent in W;. In this paper, we show that both the works 
W; and yVc satisfy the asymptotic fluctuation theorem. 
We also discuss dramatic differences between the work 
distribution functions for Wi and Wc for finite times. 

The outline of this paper is as follows. In Sec. Ullwe in- 
troduce a dragged Brownian particle model with inertia, 
and treat its dynamics using a path integral. In Sec. IIIII 
we introduce the works done in the laboratory and co- 
moving frames and calculate their distribution functions. 
In Sec. |lV]we prove the asymptotic work fluctuation the- 
orem. In Sec. |V] we discuss inertial effects in the work 
distribution functions for finite times. Finally, Sec. I VII is 
devoted to a summary and some remarks on this paper. 



II. DRAGGED BROWNIAN PARTICLE WITH 
INERTIA 

We consider a Brownian particle confined by a har- 
monic potential, which moves with a constant velocity 
through a fluid, as discussed in our previous paper [^. 
The dynamics of this particle is described by a Langevin 
equation 



cPxs 



ds 



K {Xs - VS) + (s 



(3) 



Here, we assume that is the Gaussian-white ran- 
dom force whose probability functional -Pc({Cs}) for 
{Cs}se[to,t] is given by 



PdiCs}) 



Crcxpf— / ds 
V 4a 4 



(4) 



with the normalization coefficient and the inverse 
temperature P = l/lksT), where is the Boltzmann's 
constant and T is the temperature of the heat reservoir. 
[Note that the coefficient can depend on the initial 
time to and the final time t, but such time dependences in 
C(, as well as in similar coefficients Cx and Cs introduced 
later, are suppressed.] It follows from Eq. ([4]) that the 
first two auto-correlation functions of the random force 
Cs are given by (Cs) = and (CsiCs2> = {2a/l3)S{si - S2) 
with the notation (• • • ) for an initial ensemble average. 

Now, we consider the probability functional Px{{xs}) 
for a path {xs}s£[to,t] of the particle position Xg- By 
inserting Eq. ^ into Eq. ^ and interpreting the prob- 
ability functional -P^({Cs}) for Cs as the probability func- 



tional Px{{xs}) for Xs, we obtain, apart from a normal- 
ization coefficient. 



PxiiXs}) 

= Cx exp 



; — VS m 

\ X, 

Tr a 



(5) 



with Xs = dxg/ds, Xs = d'^Xg/ds'^ and the normalization 
coefficient Cx- Here, D = ksT/a is the diffusion con- 
stant given by the Einstein relation and Tr = a/ k is the 
relaxation time in the over-damped case. For another 
derivation of Eq. ([5]) via a Fokker-Planck equation cor- 
responding to the Langevin equation, see, for example, 
Ref. [13. 

For systems whose dynamics is expressed by a second- 
order Langevin equation, like Eq. ([3]), we intro- 
duce the path integration of any functional X{{xs}) 



as /, 



(xt,xt) = {xs,pj/m) 
(xto-Xtg) = {xi,pi/m) 



Vxs X{{xs}), with respect to 

paths {xs}si^(to,t) satisfying the initial (z) condition 
{xt„,xt„) = {xi,pi/m) and the final (/) condition 
(xt,it) — {xf,pf/m). Using this notation for the func- 
tional integral, the functional average (AT ({xs}) \ over all 
possible paths {xs}s£{ta,t)i as well as averages over the 
initial and final positions and momenta of the particle is 
represented by 

r r r{xt,xt)={xf.,pf/m) 
[X{{xs})}t = dxidp, / Vxs 

J J J ixt„-xtg) = {xi,pi/m) 

^^^'^^■'^ ^({2^4) 
xPx{{Xs})f{Xi,Pi,to) (6) 

with the initial distribution function f(xi,pi,to) for the 
particle position Xi and momentum pi. The normal- 
ization condition to specify the coefficient Cx of the 
distribution functional ([5]) is given by = 1 using 

the notation ^ as well as the normalization condition 
J J dxidpi /(xijPijto) = 1 for the initial distribution 
function f{xi,pi,to). 

This finishes the introduction of our model and its dy- 
namics. In the next section IIIII we introduce the work 
done on this system and calculate its probability distri- 
bution. 



III. WORK DISTRIBUTION 



Work to Drag a Brownian Particle and its 
Distribution 



In our previous paper we considered the work W 
to move the confining potential with a velocity v in two 
frames; the laboratory frame using the particle position 
Xs and the comoving frame using the particle position 
ys = Xs — vs at time s. Based on a generalized Onsager- 
Machlup theory, we showed in Ref. [6| that the work Wi 
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done in the laboratory frame is given by ds [—k{Xs — 
vs)]v, and the work Wc done in the comoving frame is 
given by J*^ ds {-nys - mys)v with = d^ys/ds^ = 
is, leading to a difference between the work W in these 
two frames by an inertial or d'Alembert-like force —mys 
. To understand this difference in a concise way, note 
first that by the energy conservation law, the work W 
is given by the heat Q and the energy difference Ai?, 
namely by W = Q + Ai?, where the energy difference 
Ai? is the sum of the kinetic energy difference and 
the potential energy difference A[/, i.e. A_E = AC/ + 
/S.K. Here, the kinetic energy difference /S.K — /S.Kc and 
Ai^; in the comoving frame and the laboratory frame are 
given by {my1/2) - imyl/2) and - {mx^j2), 

respectively, so that wc obtain the relation 



AKr = AK, 



ds mXsV. 



(7) 



Equation ([7]) means that the kinetic energy difference 
AK depends on the frames and its frame-difference is de- 
termined by the d'Alembert-like force —mxs as a purely 
inertial effect. This frame-difference of AK also ap- 
pears in the work, and leads to the relation Wc — 
Wi — /j*^ ds mXsV. A more complete explanation for this 
frame-dependence of the work is given in Ref . [Q] , based 
on a nonequilibrium generalization of the detailed bal- 
ance condition. 

To discuss these two different kinds of work done in 
the laboratory and comoving frames simultaneously in 
this paper, we consider the work defined in general by 



ds [—k{xs — vs) — (1 — 'd)raxs 



(8) 



which gives the work W/ done in the laboratory case [d = 
1) as well as the work Wc done in the comoving case 
(i) — 0) by changing value of the parameter 7?.[2^ 
Using the functional average defined by Eq. 
the probability distribution (W) for the dimensionless 
work f3W{{xs}) is given by 



PUW,t)^iS{W-f3W{{xs}))\ 



(9) 



For later calculative convenience, we introduce a Fourier 
transformation £m{iX, t) of the work distribution function 
Pw{W,t) through the function f^(A, i) defined by 



(10) 



so that the work distribution function Pw{W) can be 
represented as 



1 



+ 00 



Pw{W,t)^^l dA£^(iA,t) e^^^. (11) 

The function (A, t) can be also regarded as a generating 
function for the work W({a;s}). By Eq. (fTO|) we obtain a 
useful identity 



£w{Q,t) = 1 



(12) 



used to determine a normalization constant later [Eq. 

mi 



B. Path Integral Analysis for Work Distribution 

To calculate the function £.^{\,t) from Eq. pU|) . we 
first note that 

r r rixt,±t)={xj,pf/rn) 
£w{\t) ^ i i dxidpi / Vxs 



j j dxfdpf f{xi,pi,to) 



X exp 



ds Ij{xg^ Xgj Xgj s) 



(13) 



by Eqs. ([5 
defined by 



and (fTO| . Here, L{xs,Xs,Xs, s) is 



1 



L[xs, is, Xg, s) — [xs 



., — vs m 

\ X. 

Tr a 



-f A/3 [k{xs — vs) + (1 — ■d)mxs] v, 

(14) 

which may be interpreted as a Lagrangian function in- 
cluding a Lagrange multiplier A due to the restriction of 
the delta function for work in Eq. © [l][2l|. Here, as 
elsewhere in this paper, the dependence of L{xs, Xs,Xs, s) 
on the parameters v, "d, etc., has not been explicitly in- 
dicated on the left-hand side of Eq. ([M]) . 

The first step to calculate the function £w{\,t) is 
to specify the most-contributing path {x*}s^^to.t] the 
path integral involved on the right-hand side of Eq. (I13|) . 
Such a special path {x*g\ s^[to,t] is introduced as the one 
satisfying the variational principle 



5 I ds L{x*,x*,xl,s) = 

'to 



with the four boundary conditions x 



to 



(15) 
■ Pi/m, 

Xf and = pf/m. In a way similar to derive 
the Euler-Lagrange equation from the minimum action 
principle in analytical mechanics ^T^, Eq. ([TS]) leads to 

d'^ dL{il,xl,xl,s) d dL{il,xl,xl,s) 



ds'^ dx* ds 

dL{x*,x*,x*, s) 



dx* 



dxt 







(16) 



for the Lagrangian function Inserting Eq. ([T] 

into Eq. we obtain a fourth-order linear differential 
equation 



4 ~* 



d^X 



1 



d^i: 



+ 



1 



rX, 







for the function x* of s, which is defined by 

xl = xl — vs + {1 — 2X)vTr, 



(17) 



(18) 
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using the inertial characteristic time r„i = m/a. 

We consider solutions of Eq. (fT7|) of the form exp(i^s). 
Inserting x* = exp(j/s) into Eq. ((T7]) we obtain the 
quadratic equation 



T^^ — I i — 2 — 



1 



1 

Tr 



= 



(19) 



for V. The solutions of Eq. (|19p are = M). , , — f_ , — ^+ 
using z/,^ defined by 



vector b|°^ can be regarded as the initial (final) posi- 
tion and velocity of the particle in the comoving frame, 
respectively. 

As the next step, we represent a path {xs\se[to,t\ as 
the sum of the most contributing path {x*}si^[to,t\ given 
by Eq. ((^ and its deviation {/!^Xs}s<^[tQ.t] defined by 



(27) 



where the variable Axg satisfies the four boundary con- 
ditions Axtg — Axj — and Aij^ = Aij = with 
Ais = dAxs/ds. Using this variable Axs, the complete 
time integral Jf^ds L{xs, is, s) of the Lagrangian 
function can be represented as 



1 

27„ 



1± ,/l-4— 



(20) 



The general solution of the fourth-order differential equa- 
tion P?)) is represented as a superposition of these special 
solutions exp(z^s) , v = , , , — v^ , namely 



Xg — Lie + 



626 - -I- 636 - 



(21) 



with constants Cj, j = 1,2,3,4. Using Eqs. ([T8|) and 
(|2ip and introducing the four dimensional vector C = 
(Ci C2 C3 6*4)^,113] we can rewrite 



xl = C^K, + vs-{l- 2X)vTr 
where the vector Ks is defined by 

/ e + \ 



(22) 



(23) 



— U, S 



The constant vector C is determined by the four bound- 
ary conditions for x* and we obtain 



where the matrix At is defined by 



(24) 



At 



/ u, tci V tn 

I e + e - 

I/, tn V in 

V, t V t 

e + e - 

I V, t ft 



-V_ to 



'.to 



-V, t 

e + 

-V, t 

-v^e + 

-u, t 

e + 



(25) 



ds L[xs, Xg, Xg, s) 



to 



= / ds L{x*,x*,x*,s) 

Jto 



J- [ ds (axs + —Axs + —Axg^ (.28) 
'^D Jto \ Tr a J 



in terms of the two variables x* and Ax*. Inserting Eq. 
([28| into Eq. ^ we obtain 



£w{\t) = Cs j J dxidp^ J J dxfdpf 



xf{xi,pi,to) exp 



ds Li^x g , X Q , X g , 



to 



(29) 



VAx, 



where Ce is defined by = /(A^o'lit )i°(o.o) 

exp[-(l/4i:>)/j*^ /j*^ds {AXs + {l/Tr)AXs + 

{raj a) AxsY'] ^-nd is independent of A. In the ex- 
pression (|29p . the contributions of the deviations Aa;^ to 
the path integral in the function £^(A,i) are included 
only in the coefficient Cs. 

Next, we calculate the quantity ft^ds L{x*,x*,xl,s) 
using Eq. ([22]), and then the function f^(A,i) given by 
Eq. (|29p. For such a calculation, using Eq. we first 
note that 



ds 



(30) 
(31) 



(z) 

and the vector B^^ is defined by 



( Xi- vto 
Pi/m - V 
X f ~ vt 
V Pf/m ~ V , 



ZVTr 



/I 



Vo 



(26) 



It may be noted that the first component Xi — vto a-nd 
the second component (pi/m) — v (the third component 
Xf — vto and the fourth component (p// m) — v ) of the 



where the 4x4 matrix 8 is defined by 



e = 




Then, using Eqs. (HH), ^ and ^ we obtain 

-x! = C^EK, -t- 2Xv 



vs 



(32) 



(33) 
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where the matrix T is introduced as 



by Eqs. ((30)) and (|3T|) . we further obtain 



r = T,„e^ + e + — X 

Tr 

\ 


o/ 



(34) 



with the relation TmV^ 



= and X the 4 X 4 



identity matrix. Using Eqs. (H!]), ([22]), jST]) and ([33]) we 

obtain 



ia/3C'^rK,KfrC 



-A(l - A)a/3w2. 



(35) 



Noting Eq. ^ and that 



ds C^QKs 



ds c^e^K 



ds 



dxZ 



v{t - to) 



to 



/' 

J to 



Xf - Xi - v(t - to), 



m 



(36) 
(37) 



ds L(x*, i*, X*, s) 



to 



1 



= aB 

4 ^ 



r(1-2A) 



Mr'' - aPXvriMf 
X{1- X)a(3v\t~to) (38) 



where the 4x4 matrix At and the vector t] are defined 

by 



(39) 
(40) 




respectively, with the 4x4 matrix $t defined by 



ds K,Kf . 



(41) 



Inserting Eq. ^ into Eq. ^ we obtain 



f^(A,i) = c'£e-^(i-^)"''"' / / dx^dp, / / dxfdpffix,,p,,to) 



X exp — 



r(1-2A) 



(1-2A) 



aPXvr}^ B 



Tr(0) 



(42) 



Equation ([42]) gives a concrete form of the function £^{X,t) for any initial distribution function f{xi,pi,to). 

The A-independent normalization coefhcient Cf in Eq. (|1^ can be determined from the condition (fT^ . and we 
obtain 



dxidpi j j dxfdpff{xi,pi,to) exp -'^aP 



B 



j(i) 



(43) 



Note that by using the condition (|T2|) we avoided to carry out explicitly the path integral included originally in the 
quantity Cg [cf. Eq. ([29]) ]. 

Inserting Eq. (j42|) into Eq. (fTTj) . and carrying out the Gaussian integral over A appearing then in Eq. (fTTj) . we 
obtain 



PUw,t) 



Cs 



y^iiraPv'^ {t-to- r2J^Af J) 



/ dxfdpf f{xt,pi,to) 



X exp < 



-a/3 



B 



{W^ - a/Jt' [t^(t - to) + {rf ~ r.J^At) bJ] } 



4a/3w2 (i - to - r2 J^At J 



(44) 



r 



where the 4-dimensional vector J is defined by 

Vo/ 



and we used the relation ni^i — 0. Equation (|44|) is 
an explicit form for the work distribution function for all 



(45) 
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time, and for any initial distribution function f{xi,pi, to)- 
Using Eq. for the coefficient Cs, the work dis- 

tribution function (|44p is properly normalized, namely 
/ dW Pni{W, t) = 1, at any time t. 

In the next two sections ITVl andrVl we discuss, using the 
work distribution function (|44|l . fluctuation properties of 
the work from the viewpoint of the asymptotic fluctua- 
tion theorem for t — > -l-oo, as well as for finite times. 



IV. ASYMPTOTIC FLUCTUATION THEOREM 



From Eq. (|47p we immediately derive 



lim ^"'(^'^) 



(48) 



for any initial distribution function f{xi,pi,to). We will 
call Eq. the asymptotic fluctuation theorem for 

work. Equation is independent of the value of the 
parameter i?, i.e. of the frame of reference (laboratory or 
comoving) or also of the contribution of the d'Alembert- 
like force to the work (0. 



The matrix At defined by Eq. ([55)1 satisfies the condi- 
tion 



lim — 

t^+oo t 



1 



to 



A* =0, 



(46) 



as shown in Appendix |^ Equation ([46]) implies that 

v{t~to) + iT]^ ~TrJ^At)Blf v{t-ta) Siudt-to- 

T^J-^AtJ t — to in Eq. ((44|) . so that the work distri- 

bution function P^{W, t) is proportional to the Gaussian 
function exp{-[W^ - aPv^{t ~ tQ)f l\\aPv'^{t - to)]} in 
the long time limit t — s- +oo, i.e. 



■ + 00 



1 



^Airafiv'^ {t - to) 

i \W ~ a(3v^ {t - tQ)Y 



X exp ■ 



AaPv^{t~ta) 



(47) 



V. INERTIAL EFFECTS FOR FINITE TIMES 



A. Slope of ln[P^{W,t)/P^{-W,t)] and the Critical 
Mass 



In contrast to the asymptotic work distribution func- 
tion (|47p , various inertial effects in the work distribution 
function appear for finite times. In this section we discuss 
such inertial effects using the function G{W, t) defined by 



G(Wt)- ^ In ^-(^'^^ 



(49) 



The function G{W, t) gives the slope of the fluctuation 
function ln[Ptu(VK, t) /Pyj{—W, t)] with respect to W, and 
satisfies 



regardless of the initial distribution function f{xi,pi, to). 
It is important to note that the work distribution func- 
tion (j47[) in the long time limit t — )■ +oo in the inertial 
case is the same as in the over-damped case. Physically, 
this is, of course, due to the finiteness of the inertial char- 
acteristic time Tm , which makes inertial effects disappear 
in the long time limit. Nevertheless, the proof of this 
equivalence is non-trivial. 



lim G{W,t) 

t — >-+oo 



1 



(50) 



by the asymptotic fluctuation theorem (|^ .[28l| 

The behavior of G{W, t) for finite times depends on the 
initial condition. To get concrete results, in this section 
we concentrate on the case of a nonequilibrium steady 
state initial condition, which can be represented by 



fix,, Pi, to) = ^W — exp<{ -/3 



[pi — mv)'^ 1 



2m 



+ ^k{x, - Vto + VTr)"^ 



(51) 



for any frame. The initial distribution function (jSip gives a Gaussian distribution for the particle initial position Xi 
and momentum pi around their nonequilibrium steady state average values vto — VTr and mv, respectively. Inserting 
Eq. ([51]) into Eq. the work distribution function is given by 



PAW,t) = ^ 
where fit is defined by 



AnafSv^ (t - to - r^J^At^) 



exp 



1 - n 



Aaf3v'^ [t - to 



^3j^[l^-a/3.^(t-to)]^} (52) 



Qt = {ri~ TrAtJ)^ [(J7 - TrAtJ) (r? - A^ J)^ + (t - to - r.^ A* j) (a^") + A^ j J {rj - A* J) . (53) 

I 



with the 4x4 matrix A^^^ defined by 

/ K 

(0) ^2 m 

- a 

V 



[See Appendix IB] for a derivation of Eg. ([5^ .] Note that 



(54) 
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the work distribution function (j52p is Gaussian with the 
average work {W) — af}v^{t — tQ) at any time because we 
chose a Gaussian nonequihbrium steady state initial con- 
dition ([?T|) . Since the work distribution function P{W, t) 
is Gaussian, G{W, t) defined by Eq. (|49)) is independent 
of W, so that we denote it by G{t)[— G{W,t)] from now 
on. Inserting Eq. into Eq. (P^j) . we obtain 



G{t) = 



1 - r^J^At J 

t-tn ^ 



(55) 



as an exphcit form of G{t). One may notice that G{t) in 
Eq. (|55p is independent of the dragging velocity v and 
the inverse temperature fi. Moreover, G{t) is positive 
for t > to because the distribution function Pyj{W,t) is 
normalizable so that the coefficient (1 — r2()/[4a/?w^(t — 
to - j'^At J)] = G{t)/['ia(3v^{t-to)] in the exponent of 
the Gaussian distribution function (j52p must be positive. 

As a first approximation to the asymptotic relaxation 
of G{t) to its final value ([SO)) , we obtain from Eq. ([SI 




G(t) 



■ + OC 



1 



FIG. 1: Time oscillation period T,n (solid line) in xl as a func- 
tion of mass m normalized by the critical mass m* = a^/(4K) 
for the mass in m/m* > 1. There is no time-oscillation for 
m/m* < 1 and the minimum of Tm is at m/m* — 2. Here, 
we used parameter values a — n = 1. We also plotted a time- 
oscillation period Tm ^ — Tm\^^Q {— 2-n 'JrnjTi) (broken line) 
for a purely harmonic oscillation in the case without the dis- 
sipation (a = 0). The time-oscillation period Tm approaches 

-(0) 



t — to — Tr 



(56) large mass limit m/m* +oo. 



meaning that the function G{t) decays to 1 by a power 
inversely proportional to the time in the long time limit 
t — » -l-oo. [See Appendix [C] for a derivation of Eq. 
(I56p .] Equation (|56p is only the first approximation for an 
asymptotic form of G{t), but already includes an impor- 
tant inertial contribution to G(t), as well as an interesting 
frame dependence of G{t). Actually, the second term on 
the right-hand side of Eq. ((56|) depends on the mass m 
via T,„ = m/a in the laboratory frame = 1, while that 
term is independent of the mass in the comoving frame 
-& — 0. Another interesting property of G(t) expressed by 
Eq. ([TO]) is that in the laboratory frame i? = 1 the second 
term on the right-hand side of Eq. ([SS]) . the t~^-decay 
term of G{t), vanishes in the case that = r™, i.e. for 
a special mass value m = o? / k. 

Perhaps the most interesting implication of Eq. ([55]) 
for G(<), although it does not appear explicitly in the 
asymptotic expression ([55]) of G(t), is the existence of 
a critical value of the mass m — m* above which G(t) 
shows a time-oscillatory behavior. In our theory, this 
time-oscillation has its origin in the time-dependence of 
X* given by Eq. (|2ip via the exponential terms exp(i'_j_t), 
etc., when the coefficient given by Eq. ([^0]) has an 
imaginary part, namely when the condition 



corresponding to a frequency uj = 
V4T,„Tr^ - l/(2r„i) \Im{v^}\, using the imagi- 

nary part Im{i/^} of v^. In Fig. [T]the time-oscillation 
period Tm (solid line) is shown as a function of the scaled 
mass m/m* for m/m* > 1. There is no time-oscillation 
of X* in the case of m/m* < 1, and the time-oscillation 
period diverges when m/m* ^1 + 0. The oscillation 
period Tm decreases rapidly as a function of mass m 
for m/m* < 2, has a minimum at m/m* = 2, and 
increases gradually for m/m* > 2. For comparison, we 
also plotted in Fig. [T] the scaled mass dependence of 
the time-oscillation period Tm '' = ^-n^m/ k (= %n\a=o) 
(broken line) for a purely harmonic oscillator with spring 
constant n. Different from the time-oscillation period 
([58)) . the period Tm^ is defined for all the masses, and 
increases monotonically as the mass increases. The 
time-oscillation period ([55]) approaches Tm^ in the large 
mass limit m/m* — s- -foo. 

It is useful to consider the critical behavior in the time- 
oscillating behavior of G(t) as due to the presence of 
two independent time scales appearing in our model: one 
characterized by Tr(= a/ k) and another by t„i(= m/a). 
These time scales r„j and are related by = 4rm* 
at the critical mass m = to*. Using these two time 
(derived from the condition ATmT/:^ - 1 > 0) is satisfied. scales, the time oscillation period T^a^ for a purely har- 



m > m = 



4k 



(57) 



We call the mass to* the critical mass in this paper, since 
a (smooth) "dynamical" phase transition takes place at 
TO — TO*. For masses m > m*, the position 5;* has a 
time-oscillation with the oscillation period Tm 



monic oscillator is given by Tm'' — 
ducing the frequencies w'"-* = 1-n /Tm'' and w, 
corresponding to the two time scales Tm^ and Tm^ re- 
spectively, the frequency lo = 211 /Tm is represented 
tj^/4 corresponding to the time- 
In this expression for the fre- 



Intro- 



Tm. — ^TTt / 1 



-1/2 



(58) 



as w = 

oscillation period 
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quency u the time oscillations occur only when the con- 
dition > uj'^/4: is satisfied. The existence of these 
two time scale and is therefore essential for the 
time-oscillatory behavior with the frequency uj, noting 
that there is no time-oscillation in the over-damped case 
containing only r^. 

In the next two subsections IV Bl and IV CI we inves- 
tigate properties of G{t) in more detail, including its 
time-oscillating behavior, for (A) the work done in the 
laboratory frame (i? = 1), and (B) the work done in the 
comoving frame (i9 = 0), separately. We will also com- 
pare those results with those for the over-damped case. 
For this purpose, we now calculate G{t) explicitly in the 
over-damped case. In our previous paper [6|, we already 
calculated the work distribution function Pw\w,t) for 
the over-damped case, which is given by 



y^AiraPv^ [t-to- t^(1 - h)] 

[W -a(3v^{t-to)y 



X exp ■ 



4al3v^t ^ to - Tril - bt)] j 
(59) 



with bt = exp[~{t~-t())/Tr] in the case of a nonequilibrium 
steady state initial distribution function f^^^Xi^to) = 
\J (in/ (27r) exp[— — vt^ + VTr)'^ /2] for the particle 
position Xi for the over-damped case at the initial time 
to-iH Using Eq. and defining, [cf. Eq. (|49l) ]. 



G^^\t) = {d/dW)ln[Pi°\w,t)/Pi"\-W,t)], we have 



Tr{l - bt) 



t-to- Tr{l - bt) ' 



(60) 



which gives G{t) for the over-damped case |21[. Note 

that Eq. §Q implies G("'(i) 1 -I- T^/(t -to- r^), 

which is consistent with Eq. (j56p , since is zero for the 
over-damped case. 



B. G{t) in the Laboratory Frame 

In this subsection we consider G{t) given by Eq. ([55]) 
— which depends on the parameter -d to specify a frame 
via At and f2t — for the work done in the laboratory 
frame, i.e. for -d = 1. In this subsection IV Bl as well as 
in the next subsection IV CI we use the parameter values 
a — K — I and set the initial time to — 0, i.e. Tr=l as a 
time unit and m/m* = Arm as the scaled mass. 

Figure [5] shows G{t) given by Eq. ([55)1 as a function 
of time t for the scaled masses m/m* — (over-damped 
case), 0.999, 2, 4, 8, 20 and 40. The graphs of G{t) aU 
converge to 1 in the long time limit t — > -|-cxi, as required 
by the asymptotic fluctuation theorem ([35]), i.e. by Eq. 

We now discuss in some detail the properties of Fig. [21 
This figure shows that G{t) for nonzero masses is always 
smaller than in the over-damped case of zero mass. In 
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FIG. 2: (color online) Graphs of G(t) = 
{d/dW)\n[P^{W,t)/P^{-W,t)\ as a function of time t 
for the work done in the laboratory frame (i? = 1) in the 
case of a nonequilibrium steady state initial condition for 
t £ [0,60]. Inset: Graphs of G{t) in a short time period for 
t £ [0,3]. Lines in these graphs correspond to parameter 
values of the scaled masses m/m* = (over-damped case), 
0.999, 2, 4, 8, 20 and 40 as indicated above this figure, and 
we used parameter values a = k = 1 (so that rr=l for a 
unit time and also m/m* = 4rm as the scaled mass, with 
m* — 1/4) and to — 0. 



the over-damped case, G{t) decreases monotonically to 
the final value 1 from -|-oo at the initial time. A similar 
behavior is still observed for small masses (e.g. see the 
graph for m/m* — 0.999 in Fig. [5] It may also be noted 
that for small nonzero masses the relaxation of G{t) to 
its final value 1 is faster than in the over-damped case 
(e.g. see the graphs for m/m* — 2 and 4 in Fig. [2]). This 
feature can be explained by the second term on the right- 
hand side of Eq. ([55)l . since the absolute value [t^ — t,„| 
of the numerator of this term is smaller for ^ — 1 than 
the corresponding over-damped value in the case of 
< m/m* < 8, using that \Tr — Tm\ < TV. Moreover, 
Fig. [21 shows that for large masses (e.g. see the graphs 
for m/m* > 4 in Fig. G{t) is smaller than 1 for long 
times, while G{t) is always larger than 1 in the over- 
damped case. This is because the second term on the 
right-hand side of Eq. ([TO)) is negative for < r,„ (i.e. 
m/m* > 4), when = 1 and t > to + Tr — t^. 

A time-oscillatory behavior of G{t) is clearly visible in 
Fig. [5] for large masses, i.e. for m » m*. To show more 
clearly the time-oscillatory behavior of G{t) for m > m* 
as opposed to for m < to*, we plotted in Fig. [3| the 
absolute value of the deviation [3011 



AG{t) = G(t) - 1 



' rr), L/ 



(61) 



t-to- Tr + Trnl9^ 

of G{t) from its asymptotic form ([55)1 as a function of 



m/m*=0.0 m/m*=0.999 m/m*= 4 - 

m/m*=0.5 m/m*=1.1 m/m*=20 - 

m/m*=0.9 m/m*=2 




FIG. 3: (color online) Linear-log plots of absolute value 
|AG(t)| of the function AG{t) = G'(f) - 1 - (r,- -r™^)/(f-r,. + 
Tm) as a function of time t for the work done in the laboratory 
frame in the case of a nonequilibrium steady state initial con- 
dition. Lines in these graphs correspond to parameter values 
of the scaled masses m/m* = 0, 0.5, 0.9, 0.999,1.1, 2, 4 and 
20 as indicated above this figure. The minima of the oscilla- 
tions of |AG(f)| for m/m* > 1 are actually zero, which is not 
indicated in this figure and Figs. |4l|6]and[7l We use the same 
parameter values a, k and to as in Fig. [2] 
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FIG. 4: (color online) Long time behavior of |AG(f)| as a 
function of time t as linear-log plots for the work done in the 
laboratory frame in the case of a nonequilibrium steady state 
initial condition. Here, we use the same parameter values a, 
K and to as in Fig. [2] Broken, dotted and dash-dotted lines 
in these graphs correspond to parameter values of the scaled 
masses m/m* = 100, 200 and 1600, respectively. Solid lines 
are fits of jAG(t)| to the function (|62|) using Table[Il together 
with the time-oscillation period (|58p . but they are visually 
indistinguishable from the graphs of |AG(t)j except for short 
times. 



time t e [0, 25] for the cases of m/m* = 0, 0.5, 0.9, 0.999, 
1.1, 2, 4 and 20. To illustrate the long time behavior of 
|AG'(i)| in more detail, we also show in Fig. [4] the abso- 
lute value |AG(t)| of AG(t) as functions of t e [0, 1000] 
for the scaled masses m/m* — 100, 200 and 1600 as 
linear-log plots. The deviation AG'(i) goes to zero when 
t —t +00 because of the asymptotic fluctuation theorem 
([50)) . In Figs. [3] and 01 it is important to note that there 
is no time-oscillation of AG(i) for < m/m* < 1, while 
we do observe time-oscillations of AG(t) for m/m* > 1, 
in agreement with a critical mass ([57)1 . above which G(t) 
oscillates in time. The decay of |AG(i)| to zero as a func- 
tion oit is faster for larger masses for < m/m* < 1 (cf. 
Fig. [3]), but slower for larger masses for m/m* > 1 (cf. 
Figs. [3] and SI). 

To check that the time oscillation period given by 
Eq. (|58p indeed appears in G(t), we fitted the data for 
AG(t) to the function 

AG{t) ae-"' sin ( + c) (62) 

with fitting parameters a, h and c in Fig. |4l The values 
of the fitting parameters a, h and c are given in Table [J 
The function ([62]) is then sufficiently close to AG(t) over 
many time-oscillation periods (except for short times) , to 
suggest that the time-oscillations of G{t) may well have 
the same origin as those in the position . Similarly for 



Fig. [21 using the fitting function ([5^ we can also check 
that the time-oscillation periods of |AG(i)| in this figure 
are given by Eq. ()58p . We fully realize that Figs. [3] and 
[4] are not enough to specify convincingly the function 
form of decay of AG(<). In Eq. ([5^ we assumed an 
exponential decay by a factor a exp(— 6t), which seems to 
fit reasonably well the data in Fig. [H However, values of 
the fitting parameters a and h shown in Table [J appear 
to vary non-negligibly if we fit data including longer time 
periods than the ones shown in Fig. [H In this sense, at 
this stage, the exponential factor in Eq. ([5^ should be 
regarded only as a convenience to check numerically the 
time oscillation period 'Tra appearing in AG(i), rather 
than claiming an asymptotic exponential decay of AG(t) 
of the form ([62| . 

C. G(t) in the Comoving Frame 

Here we consider G(^) for the work done in the comov- 
ing frame, namely the case of = 0, in which the work 
includes effects of an inertial or d'Alembert-like force. 

Figure [5] shows graphs of Gif) given by Eq. ([55[) as a 
function of time t. We chose the same masses as in Fig. 
[21 namely m/m* — (over-damped case), 0.999, 2, 4, 8, 
20 and 40 with the critical mass m* = 1/4. It is clear 
that in Fig. [digraphs of G{t) approach 1 as i ^ +oo, 
confirming the asymptotic fluctuation theorem (I48p . 
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Laboratory (1) 


100 


31.6 


-0.077 


0.021 


4.4 


Laboratory (1) 


200 


44.5 


-0.14 


0.011 


4.5 


Laboratory (1) 


1600 


125.7 


-0.32 


0.0017 


4.6 


Comoviiig (0) 


100 


31.6 


0.0038 


0.021 


4.6 


Comoviiig (0) 


200 


44.5 


0.0035 


0.011 


4.6 


Comoving (0) 


1600 


125.7 


0.0032 


0.0024 


4.7 



TABLE I: Values of the fitting parameters a, b and c for tlie function (|62l) plotted in Figs. [3] and [6] for the parameter values 
a = K = 1. 
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FIG. 5: (color online) Graphs of G(t) 
{d/dW)\n[Pn,{W,t)/Pn,{~W,t)\ as a function of time t 
for the work done in the comoving frame (i9 = 0) in the 
case of a nonequilibrium steady state initial condition for 
t G [0,60]. Inset: Graphs of G{t) in a short time period for 
t £ [0,6]. Lines in these graphs correspond to parameter 
values of the scaled masses m/m* = (over-damped case), 
0.999, 2, 4, 8, 20 and 40 as indicated above this figure and 
we use the same parameter values a, k, and to as in Fig. [2l 



FIG. 6: (color online) Linear-log plots of the absolute value 
|AG(f)| of the function AG(t) = G{t) - 1 - r,./(t - r^) as a 
function of time t for the work done in the comoving frame 
in the case of a nonequilibrium steady state initial condition. 
Lines in these graphs correspond to parameter values of the 
scaled masses m/m* — 0, 0.5, 0.9, 0.999,1.1, 2, 4 and 20 as 
indicated above this figure and we use the same parameter 
values a, k and to as in Fig. [2] 



Comparing Fig. [2]with Fig. [5l a dramatic difference in 
the behavior of G{t) in the two frames is clearly visible. 
First, a striking frame-dependence of G{t) is that for any 
nonzero mass, G{t) in the comoving frame starts from a 
finite value at the initial time io(= 0) and is always larger 
than 1, in fact going through a maximum to its final value 
1. This contrary to in the laboratory frame where G{t) 
diverges for t ^ to + and can be smaller than 1 for 
large masses and long times as discussed in Sec. IV Bl 
Another remarkable point is that, different from in the 
laboratory frame as shown in Fig. [21 G{t) converges to 
the over-damped line, much before converging to its final 
value 1, as shown in Fig. [5l This feature can be explained 
by the asymptotic form ([55)1 of G{t), whose right-hand 
side is independent of the mass m in the comoving frame 
(^? = 0), so a relaxation behavior of G{t) to its final value 



1 in this frame should be close to that of the over-damped 
case. 

Now, we discuss the time-oscillatory behavior of G{t) 
in the comoving frame. We note that in the comoving 
frame the approach of G{t) to its final value 1 is via 
oscillations around the over-damped line, contrary to in 
the laboratory frame where this approach is unrelated to 
the over-damped line. Such time-oscillations are already 
visible for large masses m » m* in Fig. [51 but to show 
them in a more magnified way, we plotted in Fig. [6] 
the absolute value |AG(i)| of the function AG{t) defined 
by Eq. as a function of time t as linear-log plots. 

Here, we plotted data for the scaled masses m/m* = 0, 

0. 5, 0.9, 0.999, 1.1, 2, 4 and 20 and for the time period 
t G [0,25]. It is shown in Fig. [S] that time-oscillations 
of AG(i) occur for m/m* > 1 but not for < m/m* < 

1. Moreover, AG{t) for < m/m* < 1 decays faster. 
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FIG. 7: (color online) Long time behavior of |AG'(t)| as a 
function of time t as linear-log plots for the work done in the 
comoving frame in the case of a nonequilibrium steady state 
initial condition. Here, we use the same parameter values a, 
K and to as in Fig. [5] Broken, dotted and dash-dotted lines 
in these graphs correspond to parameter values of the scaled 
masses m/m* = 100, 200 and 1600, respectively. Solid lines 
are fits of |AG'(f)| to the function (|62|) using Table [H and 
they are visually indistinguishable from the graphs of jAG(f)| 
except for short times. 



while AG{t) for m/m* > 1 decays slower with time, for 
increasing mass. These features are similar to those in 
the laboratory frame. 

In Fig. [7] we show linear-log plots of |AG(t)| as func- 
tions of t for longer times t S [0, 1000] and for larger 
masses m/m* = 100, 200 and 1600 than in Fig. [HI Com- 
paring this figure in the comoving frame with the corre- 
sponding Fig. [4] in the laboratory frame, we can see that 
the time-oscillation amplitudes of the function AG{t) in 
the comoving frame are much smaller than the corre- 
sponding ones in the laboratory frame, except for short 
times. This should be noted as an important frame- 
dependence in the behavior of G{t). 

The time-oscillation periods appearing in Figs. [H] and 
[7] can be checked by fitting the data again to the function 
([S^ with the time-oscillation period ([55]) . We only show 
such fitting lines for Fig. [7]using the fitting parameters a, 
b and c of Table IH Like for the fitting lines in Fig. [H the 
parameter values of a and b in Table U in the comoving 
frame also appear to vary non-negligibly for data for a 
longer time period than that shown in Fig. [T) Therefore, 
as for Fig. 31 the fitting lines in Fig. [7] should not be 
regarded as evidence for the exponential decay in the 
fitting function (|62p . However, the fits of their time- 
oscillation periods of |AG'(t)| to the function (|62|) in Fig. 
[7] are satisfactory, which suggests again that the time- 
oscillations of G{t) have the same origin as those in the 
position x*, like in the laboratory frame. 



VI. SUMMARY AND REMARKS 

As a summary of this paper, we have discussed inertial 
effects related to the particle mass m in nonequilibrium 
work distribution functions and their associated fluctu- 
ation theorems for a dragged Brownian particle model 
confined by a harmonic potential using a path integral 
approach for all times: asymptotic as well as finite. We 
considered two kinds of work: the work Wi done in the 
laboratory frame and the work Wc done in the comov- 
ing frame and we calculated the distribution functions 
Pw {W, t) for them. Using the distributions for the work 
in the different frames we analytically proved, for any 
initial condition, an asymptotic work fluctuation theo- 
rem, which has the same form in both the frames. This 
contrasts with what happens for finite times, when for 
a nonequilibrium steady state initial condition there are 
major differences between the work fluctuations in the 
laboratory and comoving frames. This was discussed, us- 
ing the quantity G(t) = {d/dW)\n[PyjiW,t)/Pyj{~W,t)], 
which approaches the value 1 in the long time limit 
t — !■ -|-oo by the asymptotic fluctuation theorem. The 
G{t) for the work Wc done in the comoving frame is larger 
than 1 at all times and converges to the corresponding 
over-damped value much before converging to its final 
value 1. On the other hand, the G{t) for the work Wi 
done in the laboratory frame can be smaller than 1 for 
sufficiently large times and masses, and the relaxation 
behavior of G{t) to its final value 1 is very different from 
that for the over-damped case, even for long times. As 
one of the significant effects for finite times, we also dis- 
cussed the existence of a critical mass m*, so that for 
the mass m > m* a time-oscillatory behavior appears in 
G{t) in both frames. 

In the remainder of this section, we make some remarks 
on the contents in the main text of this paper. 

1) We have discussed in this paper differences be- 
tween the works Wi and Wc, which originate in a frame 
dependence of the kinetic energy difference due to the 
d'Alembert-like force as we discussed in Sec. IIII Al In 
contrast to the work and the kinetic energy difference, 
the heat (as well as the potential energy difference) is 
frame-independent even in the inertial case. Note that 
the two works Wi and Wc have the same average value 
in the nonequilibrium steady state, because their differ- 
ence can be represented as a "boundary term" 

Wi - Wc = m{xt - xt„)v (63) 

depending on a difference between the two boundary val- 
ues of is at the final time s = t and the initial time 
5 = ^0 only, so that the average of this boundary term 
m(if — Xto)v is zero in the nonequilibrium steady state. 
Nevertheless, this difference m{xt-^Xtg)v between Wi and 
Wc causes dramatic differences in the work fiuctuations, 
as shown in the subsections IV Bl and IV CI of this paper. 

2) In a different nonequilibrium model described by a 
linear Langevin equation, Ref. (l6| considered the mo- 
tion of a torsion pendulum under an external torque in a 
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m 
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K 


KV 


Torsion pendulum 




I 


V 


C 





TABLE II: Correspondences between the dragged Brownian particle model described by Eq. ([3]) and the torsion pendulum 
model described by Eq. (|64p . 



fluid. The corresponding Langevin equation for the an- 
gular displacement of the pendulum at time s in this 
system is then given by 

I^ = -'^^^C9.,+Mg + <:s (64) 
ds-^ as 

where / is the total moment of inertia of the displaced 
mass, v is the viscous damping, C the elastic torsional 
stiffness of the pendulum, Mg the external torque, and Cs 
the Gaussian-white random force. For this model, Ref. 
[Tq] considered the case of a linear torque of 



Ms = 



(65) 



with a force constant /i. It is important to note that 
Eq. with the force (|65p has mathematically the 

same form as the Langevin equation ([3]) with the cor- 
respondences shown in Table [TTl Based on these corre- 
spondences between the two models, for example, there 
should be a critical value /* of the total moment of iner- 
tia above which a similar time-oscillatory behavior occurs 
in the pendulum model, like above the critical mass m* 
in the dragged Brownian particle model treated in this 
paper. 

For the pendulum system, Ref. [16j considered the 
work Wp done by the external torque Ms on the pendu- 
lum (p). This work is given there by 



ds {Ms-Mt„) 



ds 



(66) 



Using Eq. ([SS]) and the correspondences in Table [III this 
work corresponds to a quantity for our dragged Brownian 
particle model, viz. 



[ ds K,v{s — ^o)^— ■ 
Jto ds 

= Wi + Kv{t - to) 



xt - -v{t + to) 



(67) 



which is clearly different from the works Wi and Wc 
discussed in this paper. In other words, Wi, Wc and 
Wp give physically different kinds of work in nonequilib- 
rium steady states described by a mathematically iden- 
tical Langevin equation in a dynamical sense. We note 
that our Wi and Wc are consequences of the generalized 
Onsager-Machlup theory in Ref. 0]. We reserve a gen- 
eral discussion on fluctuation theorems for different kinds 
of work for a future publication. 

3) As another nonequilibrium model described by a 
linear Langevin equation, Ref. [2^ considered electric 



circuit models. In that case the system is described 
by a first-order linear Langevin equation, which has the 
same form as the over-damped Langevin equation for the 
dragged Brownian particle model. As a generalization 
of these electric circuit models, an inertial effect in the 
electric circuit can be introduced by including its self- 
induction. A generalization of the arguments of Ref. (2^ 
to the case including the self-induction, as well as a dis- 
cussion of the effects of self-induction on the nonequi- 
librium work (and heat) fluctuations, will be addressed 
in a future paper. Especially, it would be interesting 
to observe whether there is a critical value of the self- 
induction, above which similar oscillatory effects occur, 
as appear above the critical mass in the inertial case in 
this paper. 

4) The critical mass m* discussed in this paper for work 
fluctuations also appears in the dynamics of the average 
position (xs). In order to discuss this point, we note that 
taking the average of Eq. the average position (xs) 
of the particle at time s satisfies 



m- 



ds^ 



^ d{xs 
ds 



K {{Xs) - VS) 



(68) 



using {(s) — 0- Using (4 defined by Eq. (PO)) . the solution 
of Eq. (1551) is given by 



(xs) = v{s - Tr) + C"e-"+" + C'e-"-' (69) 

where the constants C and C" are determined by the 
average initial conditions {xtg) and (ito) and are given 

by 



C = - 



e + 



{{xj-v), (70) 



to 



c" = 



((ito) -«)• 



(71) 



Since the include nonzero imaginary parts for m > m* , 
a time-oscillatory behavior appears in the average posi- 
tion (xs) for masses above this critical mass m*. This 
kind of phenomenon was discussed for a damped oscilla- 
tor model , but its effect on fluctuations in a nonequi- 
librium steady state has not been discussed to the best 
of our knowledge. 
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In Ref. [6|, we discussed that in the over-damped case, 
the most probable path, which is a solution of the Euler- 
Lagrange equation for the Lagrangian function in the 
Onsager-Machlup theory, is expressed as a combination 
of forward and backward paths. This is also true in the 
inertial case, in which the most probable path is given 
by a solution of the "Euler-Lagrange" equation ^TE\i for 
A = 0. To show this, we note that the exponentially de- 
caying terms exp(— i/^ s) and exp(— 1/_ s) on the right-hand 
side of Eq. refer to a forward path. We can also 

introduce the corresponding backward path, as a com- 
bination of exponentially divergent terms exp(i>^s) and 
exp(^'_s). A combination of these forward and backward 
paths gives then the most probable path {xl}selto,t] for 
A 0, i.e. Eq. 

5) There is still the open question of an analytical dis- 
cussion of the asymptotic form of AG{t) with the time- 
oscillations shown in Figs. El SI and [T] In this paper 
we only analyzed AG{t) numerically by fitting it to the 
function (|62p . but in principle, such analytical informa- 
tion on AG{t) is contained in the general form (|55)) of 
G{t). 

6) We have considered the asymptotic fluctuation the- 
orem for work in this paper. We now address very briefly 
its connection with other fluctuation theorems. 

(6a) One of the other fluctuation theorems is the tran- 
sient fluctuation theorem j2i]- This fluctuation theorem 
was already derived and discussed for a dragged Brow- 
nian particle model with inertia in Ref. Q. There, we 
derived transient fluctuation theorems, not only for the 
same works as those in this paper, but also for an energy 
loss by friction. Different from the work, the distribution 
function for the energy loss by friction does not satisfy 
an asymptotic fluctuation theorem. 

(6b) Another important fluctuation theorem is the ex- 
tended heat fluctuation theorem In Ref. @ 



we gave a simple derivation of this fluctuation theorem, 
based on the assumptions that (A) a correlation between 
the work and the energy difference at time t (as well as 
a correlation between the energies at the initial time to 
and the flnal time t) disappears in the long time limit 
t +00, (B) the work satisfies the asymptotic fluctu- 
ation theorem, (C) the work distribution function ap- 
proaches a Gaussian distribution asymptotically in time, 
and (D) the distribution function Pe{E) for energy E 
is canonical-like, namely Pe{E) « exp(— /3ii^) for E > Q. 
The same derivation could be applied to all models which 
satisfy these four conditions (A), (B), (C) and (D). In 
particular, using this derivation, one can derive an an- 
alytical expression for the asymptotic heat distribution 
function itself, as well as the extended heat fluctuation 
theorem not only for the over-damped case, as was done 
in Ref. 6], but also for the inertial case. 
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APPENDIX A: ASYMPTOTIC PROPERTY FOR 
THE MATRIX At 

In this Appendix we prove Eq. ([46)) for the matrix At. 
To show this equation in a simple way, without losing 
generality we take the origin of time at (io + i)/2 so that 
the initial time is given by to = only in this Appendix. 

To consider the structure of the matrix At defined by 
Eq. ([55)1 in the long time limit t ^ +cx3, we first calculate 
the asymptotic form of the matrix F^tF, which is an 
essential element of the matrix At. For this purpose we 
note 



)s 



2u s 



1 

-2i/ s 



V 



1 e' 



Inserting Eq. (jA.ip into Eq. ([41]) and using the relation = ~i we obtain 

/ sinh{2;-'_^ i) sinh[(i>^ )t\ sinh[(^'_^ )i] 
o7] 



2v^ v^+v_ 
sinh[(iA +i/_ )i] sinh(2i/_t) 



t 

sinh[(i/^ — z/_ )t\ 



sinh[(iA — )t] 



2v_ 
t 



h(2i/_t) siiih[(iA +!/_ )t] 



\ 



t 



sinh[(^'_|_ )t] sinh[(i^_^ )^] sinh(2;-'^ t) 



J 



with the hyperbolic function sinh(a;) = [exp(x) — exp(— a;)]/2. Equations ([M]) and (IA.2P lead to 

"ft O2 



F$tF = 



O2 O2 



(A.l) 



(A.2) 



(A.3) 
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where O2 is the 2x2 null matrix, and the 2x2 matrix '^t is given by 



smh{2v_^t) TT+IT sinh[(j/_^ + v_ )t] 

^ sinh[(j/^ + iy_)t] V sm\i{2v_t) 

2u, t 2v V_ (1/ +iy )t 

Iv J^v )t 2v t 

6"+* \ / ^4 U^T^ \ /e"+' 



(A.4) 



Here, we used the positivity Re{vj^} > of the real 
part of 14 (assuming a nonzero mass m 7^ and a 

nonzero spring constant k 7^ 0) and also sinh(af) '"^^i"" 
(1/2) exp(at) for any number a with the positive real part 
Re{a) > 0. 

Second, we obtain a simplified form of the matrix A^T^ 
in the long time limit, which is another essential element 
of the matrix Aj. Noting again that the real part of the 
number Vj^ is (strictly non-zero) positive and the initial 
time is given by = —t, we obtain the asymptotic form 
of the matrix At defined by Eq. ([25|l as 



Here, using Eqs. (|A.4[) and (|A.10[) the non- vanishing 
matrix elements of the matrix (|A.11[) is given by 



A 



(2)-l 



(2)_i 



where we us ed + 1'_ — a/m and 
(jA.lip and (|A.12p we obtain 




At 



for the long time limit t 



(1) 



O2 



(A.5) 



null matrix, and A. 



U) 



J 



A. 



(1) 



A. 



(2) 



+00. Here, O2 is the 2x2 
1 , 2 are defined by 



(A.6) 
(A.7) 







u, t 

e + 










( 








u^e"-' ) 








lim At = — I „ „ „ 
a [ K 

m 



(A.13) 



Equation (jA.lSp shows that the matrix At approaches a 
time-independent constant matrix in the long time limit 
t +00. Therefore, the matrix Kt/{t — to) approaches 
the 4x4 null matrix in the long time limit t +00, 
implying that the condition (|^ is satisfied. 



From Eq. (jA.SP we derive 



-1 t— >+oo 



A 



(2)-l 



A 



(1)-1 



where A^^^ ^ and Ap-* ^ are given by 



1 



A. 



(2)-l 



1 



V, 



-1/ t 



~v, t 

V e + 



-e + 



(A.8) 



. (A.9) 
(A.IO) 



APPENDIX B: WORK DISTRIBUTION FOR 
THE NONEQUILIBRIUM STEADY STATE 

In this Appendix we give a derivation of Eq. (|52p for 
the work distribution function P(VK, t) in the case of the 
nonequilibrium steady state initial condition ([5T|) . 

First, we note that the initial distribution function ((5T|) 
can be written in the form 



/(xj,pi,io) 



27r 



exp 



B 



A(o)B« 



(B.l) 



Eq. (|A.8[) give an asymptotic form for the matrix A^ ^ . 

Finally, using Eqs. dM]), (|A.3p and (|A.8p we obtain the 
asymptotic form of the matrix At as 



At 



A 



(2)-l 



using Eqs. ([26)1 and ([541) . Equation ()B.ip means that 
the initial distribution function f{xi,pi,tQ) is Gaussian 

for the components of the vector B^^-*. Using Eq. (IB.l 
and again the vector B^j 



if 

given by Eq. ([^5)) . the work 



(A. 11) distribution function ([^^ can be represented by 
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Km 



X exp 



B 



(1) 



a[3v {q^ - TrJ^At) -W + al3v^{t - to) 



4a/3w2 (t-to- rSJ^At J) 



> (B.2) 



where we used the relation dxidpidx fdpf = m^dB^^-' because of Eq. 
Now, we note that 



26a'^x + 6^ 



X + & ( aa^ +C] a 



aa^ +C 



X + 6 ( aa^ +C] a 



1 - a^ (^aa^ + a 



(B.3) 



for any n-dimensional vector a and x, any scalar b, and any n x n symmetric matrix C existing the inverse matrix of 
aa"^ + C. Applying Eq. (|B.3p to Eq. (|B.2p for the case of x = B|j.\ a = al3v{r] - t^AjJ), b = al3v^{t - to) - W and 
C = {al3vY{t - to - r^j'^AtJ)(A(°) + At), and carrying out the integral over B-l''(= x) in Eq. (IB.2I1 . we obtain 



Km 



Ceml3 

4Trv V 7ra/3 [t - to - T^3^At3) 



nil 



d5c exp ■ 



(a^x + 6)2 + x^Cx 
'4a/3t)2 - to - r2 jTAt J) 




1 - a^' (^aa^ a 



4a/3w2 (t - io - T2J2^AtJ) 

l-»t 

'4al3v^ {t-to-T^J^AtJ) 



[W -af]v^{t~to)]' 



(B.4) 



using a normalization constant C^, and Eq. (|53p for fit. 
The constant C^, in Eq. (|B.4p can be determined by the 
normalization condition J dW^ P{W,t) = 1, which and 
Eq. dEll) yield Eq. dH]). 



APPENDIX C: ASYMPTOTIC FORM OF G{t) 

In this Appendix we give an argument to derive Eq. 
^ for G{t). 

The essential point to derive Eq. (jSG]) for G{t) is the 
asymptotic form (|A.13p . or equivalently 



lim At = 2 

t — >-\-OG 















1/Tr 











(C.l) 



using the relations t,. = a/ k and t„i — m/a. Using Eqs. 
pl)l . (|i5)) . (|5i)) and (jCll) we obtain 



lim J^AtJ 2/rr, 



lim [ri - r,.At J) = 

i— *+oo 




l/r^ 



(C.2) 
(C.3) 



^limJA(")+At) = 2| I I |(.C.4) 

T„ 

Equations ((U?2l) . (IU?3)) and ([Till) lead to 
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2{t -to- 2Tr)Et 



-to- 






-A.)l 


-1 


t — Tr 


-'d 


-Tr 




--d 


TrTm 




£ 




-Tr 


-1} 


Et — Tr 






£_ 









(C.5) 



with 2f = 2 - to - + ^?V,„). By Eqs. jUj]) From Eqs. JMl), ([Cj) and ([ClI we derive Eq. JM]) 

and (|C.5p . we obtain 
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